Constraints on Topological Defect Formation in 
First-order Superconducting Phase Transitions 
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In this work we address the impact of a cubic term addition to the Ginzburg-Landau mean- 
field potential, and study the consequences on the description of first order phase transitions in 
superconductors. Constraints are obtained from experiment and used to assess consequences on 
topological defect creation. No fundamental changes in either the Kibble-Zurek or Hindmarsh- 
Rajantie predictions are found. 



I. INTRODUCTION 

The following work is based on the research reported in 
Ref. 01 . It pursues the objective of empathizing the anal- 
ogy between accessible condensed matter systems and the 
currently accepted framework for the evolution of the 
Universe, notably the symmetry-breaking phase transi- 
tions 0)01 has undergone after the Big-Bang. A key as- 
pect to this comparison is the creation of topological de- 
fects, frustrations of the unbroken phase within the bro- 
ken one, arising from the continuity of the order param- 
eter values. These are generally categorized according 
to the homotopy group of the quocient of the unbroken 
symmetry groups to the broken one and which enables 
for comparison of different physical phenomena. These 
objects can appear as magnetic monopoles or point-like 
defects, cosmic strings, vortices or flux tubes, magnetic 
domain walls or textures. Besides its mere aesthetical 
value, this analogy can provide a powerful probe into the 
early stages of the evolution of the Universe, since direct, 
hands-on experimental tests are unattainable: the exis- 
tence of more accessible systems that exhibit a formally 
similar behavior could provide crucial clues to many cos- 
mologically relevant issues. 

These "cosmology in the laboratory" experiments can 
be found in various systems, ranging from vortices in 
superfluid phase transitions of *He and '^He (see e.g. 
Ref. 0, 0]), which exhibit common features with cos- 
mic strings |y| , to liquid crystals undergoing an isotropic- 
nematic phase transition 0, 0. Polymer chains were 
shown to also possess analogous thermodynamic and 
transitional behavior 0. However, most of these sys- 
tems lack the existence of a quantity analogous to the 
magnetic field, which could be a key player in the early 
evolution of the Universe and formation of structure. For 
that reason, superconductors are a case of special in- 
terest. These comprise phase transitions involving a lo- 
cal gauge symmetry-breaking process, during which the 
photon acquires a "mass" and, therefore, a penetration 
length, giving rise to the Meissner effect: the expulsion of 
the magnetic field from a superconducting material, with 
formation of shielding "supercurrents" on its surface. 

This symmetry breaking originates topological defects 



which are known as flux tubes or vortices, lines of non- 
null magnetic field trapped inside the superconductor. 
Experiments targeted at observing defect densities in 
high-Tc materials [Tol | were not in accordance with the 
density predictions of the Kibble-Zurek (K-Z) mecha- 
nism 21. This, however, is to be expected, since the 
K-Z prediction should be accurate only for global gauge 
symmetry breaking, when the geodesic rule for phase an- 
gle summation is valid. A new defect generation mech- 
anism, based on a local gauge treatment by Hindmarsh 
and Rajantie (H-R) leads to an (additive) predic- 
tion. This, although well below the first Carmi-Polturak 
experimental sensitivity, is in reasonable agreement with 
the second. 

As a starting point for this research, we note that the 
above experiments were both conducted in type-II ma- 
terials, which exhibit a higher critical temperature and 
are therefore easier to manipulate, leading the current 
trend in experimental superconductivity. These materi- 
als display a second order phase transition, with no re- 
lease of latent heat. On the other hand, Type-I materials 
are metastable, showing different responses to a magnetic 
field when in normal-superconductor or superconductor- 
normal phase transitions. In this work, we try to account 
for this more elaborate behavior and to estimate to which 
extent it affects the defect density predictions for type-I 
superconductors. 

Type-I and type-II superconductors are commonly dis- 
tinguished according to their Ginzburg-Landau (G-L) pa- 
rameter K = A/^, the ratio between the magnetic field 
penetration length A and the coherence length ^ of the 
order parameter. In the presence of a gauge field A, these 
characteristic length scales are obtained from the free en- 
ergy density 
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+ V{^) + (V X A) 



(1) 

where fl is the sample's magnetic moment, rUe is the elec- 
tron mass and is the order parameter. The G-L po- 
tential is usually written as 0| 
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FIG. 1: Characteristic Potential curves. 



(2) 



where a is assumed to be linear with temperature, a = 
a'{t— 1), t = T/Tc, a' and (3 are constants, and Tc is the 
critical temperature. One obtains 



(3) 



In the normal-to-superconductor phase transition, a 
term similar to — 7|<i>p also arises if one takes into ac- 
count gauge field fluctuations producing 
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(6) 



This result enables a first order phase transition for all 
values of k. Thermal fluctuations Q| and non local BCS 
effects ^3 describe crossover behavior between first and 
second order transitions. 

In the following, we adopt a potential of the form of Eq. 
and constrain 7(T) based on experimental data. The 
results are then compared with both TFT 1-loop radia- 
tive corrections and the the results of Ref . ^3 (valid only 
at temperature close to Tc, that is, t 1). The intro- 
duction of 7(r) in the potential (0) can produce changes 
in both the K-Z and H-R defect generating mechanisms. 
Also, a possible nucleation suppression due to the po- 
tential barrier can significantly reduce the number of ob- 
served defects. The obtained constraints on 7(T) are 
used to access the impact on these claims. 



II. TEMPERATURE SENSITIVITY BOUND 



and 



(4) 



The coherence length at zero temperature is — 
h/ \/2meOi' , with k ~ ^/p. The transition is second order 
if K > l/-\/2, and i^o is typically less than ^ 0.04 /im; 
for K < 1/ v^, the transition is first order and typi- 
cally greater than 0.08 ^m. It is commonly accepted 
that, if there is no applied magnetic field, one always 
has a second-order phase transition, for all values of k. 
First order transitions arise from the external field term 
in Eq. if the sample has a characteristic dimension 
I > X. This degeneracy of the phase transition a,t H = 
is, however, only verifiable to the current experimental 
sensitivity, and it can be argued that there is some yet 
undetected intrinsic metastability, regardless of the ap- 
plied magnetic field. 

Bearing in mind the analogy between condensed mat- 
ter and cosmology, wc now briefly look at phase tran- 
sitions in high energy physics. In thermal field theory 
(TFT) a first order phase transition arises due to 1-loop 
radiative corrections to a potential similar to that of Eq. 

a barrier between minima of the potential is created, 
as these corrections give rise to a cubic scalar field term 



V{^) = -7|$|' 
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(5) 



where "f{T) = (\/2/47r)e^T 12]. As before, /3 is assumed 
to be constant and a — a' (t — I) to depend linearly with 
temperature. 



A superconductor undergoing a first-order phase tran- 
sition crosses different supercritical fields, displaying a 
metastable behavior, as shown in the phase diagram 
of Figure ^ The superheating curve is given by the 
condition |^ = ^ = 0, for $ 7^ 0, equivalent to 
a = 97^/16/3. The supercooling curve is given by the con- 
for $ = 0, corresponding to a = 0. The 
(unobservable) critical curve is given by V{0) — V{^c) 
and '^^^^"^ = 0, where <I>c is the non- vanishing minimum 
of the potential. This corresponds to a = 7^/2/3. 

Assuming a = a'{t — 1) and "f{t) — S t, we obtain for 
the superheating curve 
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and 
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4 a'/3 



16^ 



(8) 



Due to the presence of the cubic term in the potential of 
Eq. (O, the superheating curve shows a zero-field shift 
in temperature from Tc by {96'^ /16a' (3)Tc. This shift, 
if detected, would indicate an intrinsic metastability, in 
the sense that it does not depend on the existence of 
an applied field. Since such temperature shift has not 
yet been signaled, the current experimental temperature 
sensitivity being Atcxp 10~^ [ig, a bound on the slope 
of 7 is 
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TABLE L Critical properties of Sn and Al 



Material 


Tc (K) 


Hc{Q) (G) 


Co (m?") 


A (nm) 


Sn 


3.7 


309 


0.23 


34 


Al 


1.2 


105 


1.6 


16 



TABLE IL Derived quantities and bounds for 5 



Material 


Sn 


Al 


a' (J) 


1.15 X 10"^^ 


2.38 X 10"^^ 


P (J.m^) 


4.72 X 10"^'' 


2.16 X 10~^® 


a' (eV^) 


3.61 X 10"^ 


7.45 X 10"^ 


P 


9.45 X 10"* 


4.32 X lO"'' 



Ail 

16 a' (3 ^ 



(9) 



The supercooling transition still occurs a,t t — 1, the 
critical temperature; this is natural, since it is determined 
solely by a = (neglecting a smaller order correction to 
a, M)- 



bound 


tsh shift 


S{eV) 


tsh shift 


S{eV) 


hah 


0.25 


1.23 X 10"^ 


0.25 


1.20 X 10"" 




10~^ 


7.78 X 10"'' 


10"^ 


7.57 X 10"'' 


Ref. fl3] 


5.19 X 10"^ 


1.77 X 10"^ 


2.32 X 10"" 


3.64 X 10"'^ 


TFT 


2.92 X 10"^ 


1.33 X lO"** 


3.24 X lO"** 


4.31 X lO"'^ 



III. SUPERHEATING PERTURBATION 
BOUND 

To obtain further constraints, one looks more carefully 
at the mechanism relating the presence of an applied field 
with metastability. According to Ref. ^3li this arises 
due to the contribution of the magnetic moment to the 
free energy and hence it depends not only on the 
material (that is, on k), but also on its shape and dimen- 
sions. The supercritical fields and the value of k have 
commonly been obtained from experiments with micro- 
spheres of type-I materials. Table |2 indicates the critical 
properties of two of these, Sn and Al. To assess the influ- 
ence of the cubic term of Eq. © , we reproduce Ref. [T^ 
calculations, including the cubic term in the potential. 
For a small sphere of radius a, the magnetic moment is 
given by [l3 



h1,^{t) = ^L=h^Mi^~t') , (12) 

where /i|^(0) = Hc{0)/ Hc{t), an approximation valid only 
close to Tc 24). 

The calculation of the supercooling field implies the 
evaluation of a second derivative at the origin, and clearly 
the presence of a cubic term in the potential has no effect, 

^(7G<i>') = 67G<i>^0 . (13) 

Since no shift of the superheating curve due to 7 has 
been detected, we must have 7g ^ 1- This constraint 
fails for relative temperatures in the range 
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3A ^a$o 
cotn 



3A^ 



H 
8^ 



(10) 



where $0 = ^/^oc and <1>^ = rUeC^ /iire'^X'^ . 

After a computation presented in the reduced su- 
perheating field hsh = Hsh/Hc is given by 



as 



IG 



(11) 



where 70 = 37/2-\/|a|/3 is defined to be dimensionless, 
while /i"^ is the unperturbed superheating field, corre- 
sponding to 7 = 0. 

Measurements [n m m m m were obtained with 

colloidal dispersions, and the statistical error due to the 
size distributions of the microspheres do not allow for 
a direct fit of 7g from hsh{t) data. The measurement 
reported in Ref. [23| used single microspheres, but non- 
locality and impurities did lead to a large theoretical un- 
certainty in 



9 52 

1 < t < 1 

4a'/3 

This interval is vanishingly small if 



(14) 



9^ 

4a'/3 



< 1 



(15) 



which is a weaker bound than the one of Eq. @ . 

Al and Sn show maximum critical fields of order 10^ G, 
and the shift between hsh and is smaller than 10"^ G. 
Since this is well below the sensitivity of measurements 
20, 21, 22, 23], we consider only Eq. © and 
drop the bound of Eq. (fT^ . Notice that, even if the 
"exclusion" interval" (|14|l is non-negligible, no "spikes" 
should appear in that region of the H — T diagram, due 
to the smallness of the field values there. 

Table provides a comparison of the bounds on 5 
with the cubic term arising from 1-loop corrections in 
TFT and the prediction of Ref. This is achieved by 
computing the slope of ^{t) from 7(T) = (•\/2/47r)e'^T, 
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obtaining 6 = (\/2/47r)e^Tc. Notice that 7(t) as a func- 
tion of the reduced temperature is material dependent, 
although 7(r) is not. The quantities a' and f3 are also 
included. 

We notice that unit conversion is not direct, but 
achieved through a multiplicative factor mg! since the 
dimension of the scalar field in G-L theory is [$^] = L~^, 
its square representing a density, while in field theory 
[$] = L~^, the dimensionality of 7 depends of the theory 
at hand. For comparison sake, we have chosen [7] = L~^. 
Hence, the definition of 7 is changed with respect to 
the free energy potential of Eq. through a conve- 

nient TTie factor; the electron mass determines the con- 
version as it is absent from the kinetic term of the La- 
grangean density of field theory, 9^<I>(9'^$, but present 
in the corresponding condensed matter free energy term, 
(?i^/2me)V^<i>; equivalently, one can loook at the coher- 
ence length: S^pj, = 1/a' vs. = Ti^ /2mea' . 

According to Ref. 0, in the absence of an apphed 
magnetic field, momentum fluctuations of the gauge field 
have an expectation value derived from the equipartition 
theorem. Integrating over momentum space (with a cut- 
off A of the order oi S^^^), one gets 



IT 



8Mo^%/7i7*oTc|$| 
Tic 



(16) 



topological defect per area is assumed. Since the char- 
acteristic scales ^^^d A are obtained by linearizing the 
G-L equations close to Tc, thus neglecting the 7-cubic 
and /3-quartic terms contribution, we expect no signifi- 
cant changes to ^0 and, hence, to this prediction. 

The defect density predicted by the H-R mecha- 
nism for a thin slab of width Lz is given by n ~ 
(e/27r)T^/^LJ^''^^~^, where ^ 27r/fc is the domain size 
after the transition and is related to the highest wave 
number k to fall out of equilibrium. This is obtained 
from the adiabaticity relation 



duj{k) 



dt 



(18) 



For an underdamped dispersion relation, one has 
Lu(k) = ^ -\- mi and the photon mass is given by 



= 2e^|<i>|^ = ~2e^a//3 . This leads to a defect pre- 

— 1/3 

diction n oc Tq . Recall that the penetration length 
A expresses a non-null photon "mass" arising from the 
spontaneous symmetry breaking that occurs during the 
normal-to-superconductor phase transition. The pres- 
ence of a cubic term in the potential will change the 
photon mass, as the true vacuum of the broken phase 
shifts to 



Hence, from this result arises, asides from an unim- 
portant correction to the scalar field mass (or coher- 
ence length), a more relevant (negative) cubic term, 
—8^o{e/fic)^n^oTc\^\^ . At zero field, this produces a 
shift in the superheating temperature of 



At = 7.25 X IQ-^^T^H^iOf^^ 



(17) 



with Hc{0) expressed in Gauss and ^0 in M*^- For Sn, it 
requires a temperature sensitivity of 10~^K, well below 
current possibilities. Al, however, requires just a sensi- 
tivity of 10~^K, attainable if one employs state of the 
art relative temperature measurement techniques. 

One can see that, for both materials, the slopes of 7 
predicted by TFT and Ref. have similar magnitudes 
^ 10~'^eV. This is a confirmation of the underlying anal- 
ogous mechanisms: one can view the thermal averaging 
of the gauge field in condensed matter as equivalent to 
finite temperature vacuum polarization, given by renor- 
malization of 1-loop Fcynmann diagrams. 



IV. TOPOLOGICAL DEFECT FORMATION 

We now discuss possible implications of the presence 
of this cubic term in the mean-field potential: the K-Z 
mechanism predicts a topological defects (vortices) den- 
sity of n ~ S.o'^iTo/TqY, where tq = nh/16kBTc is the 
characteristic time scale, given by the Gorkov equation, 
Tq is the quench time, and v is a critical exponent. One 



-37+ \/-16a/3-t-97^ 



(19) 



However, since 70 = ^1 /'2\/\a\f3 <C 1, the effect is, as in 
the K-Z case, too small to affect the H-R prediction in a 
significant way. 

There is also a non-vanishing probability of the order 
parameter to quantum tunnel from the false vacuum to- 
wards the true one. The rate of transition per unit vol- 
ume and time is given, in the thin wall approximation, 

mm 



VAt 



5*3 



2tt T 



3/2 



'S3/T 



2tt 1 79(T) 



(20) 



(21) 



is the Euclidean action, and e(T) the "depth" of the true 
vacuum. For the thin wall approximation to be valid, 
one assumes that the barrier's height is much greater 
than this "depth" e, which implies that 7 is comparable 
to a and (3, that is, 7g ^ 1- Since, as shown above, the 
current temperature sensitivity of 10"^^" only allows for 
7g < 10~^, the approximation breaks: the field should 
always tunnel through the potential barrier, that is, with 
a probability close to unity. Therefore, topological defect 
production is unsuppressed, and one does not need to 
be concerned that these may not have sufficient time to 
nucleate within the measurement's resolution time. 
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V. CONCLUSIONS 

In this work we have examined a description of a type-I 
superconductive phase transition including a cubic term 
in the G-L mean-field potential, inspired both by anal- 
ogy with TFT 1-loop corrections and gauge field thermal 
averaging 13]. Superheating field and temperature con- 
straints impose bounds on this cubic term so that its 
effects are small compared to that of other parameters in 
the G-L potential. In particular, there is negligible im- 
pact on the defect density predictions of K-Z 2, 6] or H-R 
with no suppression or slowing down of defect pro- 
duction due to nucleation suppression arising from the 
induced potential barrier between false and true vacua. 

In the absence of any sizeable perturbations, the H-R 
prediction for topological defect density in type-I super- 
conductors should be reduced by a factor of 10 — 100. 
However, it has been suggested that defects nucleated 
in type-I materials survive significantly longer than in 
type-II |23|: for type-I superconductors, their lifetime is 
expected to be of order 10~^ seconds. This increases the 
possibility of measuring the created defects before they 
disappear, possibly compensating for an inferior net num- 
ber. 

Finally, we point out that the superheating tempera- 



ture shift induced by a cubic term, derived either from 
Ref. pSj or from TFT 1-loop corrections, increases with 
decreasing n {Msh{Sn) - 10"^; AU.(A/) - 10"^). 
Thus, possible searches for a cubic term should be con- 
ducted with extreme (k ^ 1) type-I materials like a- 
tungsten, which displays Tc — 15.4 ± 0.5 mK , He = 
1.15 ± 0.03 G. Also, much more precise bounds could 
be attained by using a DC SQUID to measure the shift 
in the supercritical field, since this device currently pos- 
sesses a sensitivity of 10~^cI)q/\/Hz, or 10~^G over a 10 
/im grain diameter. 
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